Nuclear dynamics in time-dependent picture 
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Abstract. Using the time-dependent theory of quantum mechanics, we investigate nuclear electric 
dipole responses. The time evolution of a wave function is explicitly calculated in the coordinate- 
space representation. The particle continuum is treated with the absorbing boundary condition. 
Calculated time-dependent quantities are transformed into those of familiar energy representation. 
We apply the method to a three-body model for ^'Li and to the mean-field model for ^^O, then 
discuss properties of £1 response. 



TIME-DEPENDENT METHOD FOR NUCLEAR RESPONSES 

Atomic nuclei exhibit a variety of responses to different experimental probes; Coulomb 
and nuclear excitations, spin- and isospin-dependent probes, high- and low-energy reac- 
tions. In order to investigate properties of these nuclear excitations, it is useful to study 
response function for a probe (external) operator of interest. Normally, these quantum 
many-body problems are theoretically studied using the energy representation, in which 
calculations are carried out by either diagonalizing the Hamiltonian matrices or solv- 
ing differential equations with a fixed energy. These approaches are advantageous when 
we are interested in states within a limited energy range. However, when we calculate 
response functions in a wide energy range, the energy representation may not be the 
best choice. Occasionally, the time representation provides more efficient numerical ap- 
proaches to quantum problems. In this paper, instead of taking the time-independent 
(energy-dependent) framework, we present time-dependent approaches to nuclear re- 
sponse functions, and show their advantages and usefulness. 



Time-dependent equation with absorbing boundary 

The quantum mechanical problems are often described by the eigenvalue equation, 

H\'i') =£{'¥). (1) 
Since this stationary Schrodinger equation is derived from the time-dependent equation 

ih^^\^'{t))=H\^'{t)), (2) 

these two equations should provide the same information if we could solve Eq. ([T]) for 
all the energies E and solve Eq. © for an infinitely long period of time. However, in 



practice, we need to truncate either the energy range or the period of time propagation. 
Therefore, these two approaches have advantages and disadvantages, complimentary 
to each other. The stationary equation, Eq. ([T]), is suitable for precise calculation of 
eigenstates in a small energy range. In contrast, the time-dependent one, Eq. 0, will be 
superior in practice, when we are interested in a bulk property of response function 
for an extended energy range. The energy resolution obtained in the calculation is, 
roughly speaking, inversely proportional to the period of time propagation; /S.E ^h/T.ln 
general, longer the time propagation is, more precise energy the calculation can predict. 
However, for resonance states with positive energies, the time period T can be safely 
truncated without much loss of the accuracy, because the states possess a finite width. 

For the positive-energy states, we must deal with a problem of the boundary condi- 
tion. For the time-dependent framework, since the asymptotic form of the wave func- 
tion becomes a superposition of a variety of momentum eigenstates, solving the time- 
dependent Schrodinger equation © with a proper continuum boundary condition seems 
to be a very difficult task. However, the use of the absorbing complex potential provides 
a practical approach to the problem ||lll2.B0]- We use the complex potential to im- 
pose an approximate outgoing boundary condition on calculated time-dependent wave 
functions. This is simply done by replace Eq. Q by 

m^^l^'it)) = {H -iW{r))\^'{t)), (3) 

where —iW{r) is a coordinate-dependent imaginary potential. Suppose we calculate the 
wave function \^{t)) as a function of time t with an initial state |^in) that is localized in 
space. W{r) must be zero in the region of ^in(?^) 7^ 0. while it is either zero or positive 
outside of this region. In addition, VM^(r) should be small to prevent the reflection, 
but simultaneously W{r) should be large enough to absorb all the outgoing waves. As 
long as these conditions are satisfied, in the region of ^in(?^) 7^ 0, the wave function, 
^(r, r), obtained by solving Eq. ^ is identical to the solution of Eq. Q with the 
outgoing boundary condition. Since small reflection of the waves by the potential —iW 
is inevitable, this boundary condition is approximate, however, we have confirmed that 
the approximation is good enough to produce results indistinguishable from those of the 
exact boundary condition [i^i^CD- 

APPLICATIONS 
Three-body model: Coulomb breakup of ^^Li 

In this section, we apply the time-dependent method to Coulomb breakup reaction 
of the two-neutron-halo nucleus, '^Li. In the perturbative regime, properties of the 
Coulomb excitation are primarily determined by dB{El)/dE values. Strong soft E\ 
excitation in ^^Li were previously reported at MSU |8], at RIKEN 11], and at GSI 
However, these data seem to be inconsistent with each other. Recently, the experiment 
has been performed at RIKEN with much higher statistics I1II1I . The new data indicate a 
strong El peak around Ey^ = 0.6 MeV. 



The B{E1) strength distribution is defined by 



= 1 1 5 (£ - £') I i^E' \^o) I ' , (4) 

where |^o) and l^^) are the ground and excited states of energy E, respectively. is 
the £1 operator with a recoil charge. This quantity is rewritten in the following way: 
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Thus, B(El) strength distribution can be calculated using the time propagation of the 
state I'i'^it)). The initial state is given by |^^(0)) = Mi^|^o) in Eq. (El). It should 
be noted that, in the time-dependent method, we do not need to construct the energy 
eigenstates, \^e)- 

We adopt the following Hamiltonian for ^^Li, neglecting the recoil of the ^Li core. 

2 2 

^ " " ^'^2 + ync{n)+Vncir2)+Vnn{\ri-r2\), (7) 

where m indicates the neutron mass. The interaction potential between neutron-^Li (Vnc) 
and that of neutron-neutron (V,„j) are both in the Woods-Saxon form with radius of 2.3 
fm and the diffuseness of 0.6 fm. Occupied orbitals, ^^(ri) and ^sifi), are excluded from 
the model space (H PHP), using a projection operator, 

P^{1- - |(^),(2))(<|),(2)|}. (8) 

The ground state is assumed to be a spin-singlet state, thus must be symmetric in the 
coordinate space, ^0(^,^2) = ^0(^25^)- The depth of Vnc potentials are adjusted so as 
to reproduce the two-neutron separation energy of 0.3 MeV. Numerical calculations are 
performed using the partial-wave expansion with respect to the n-C coordinates, ri and 
r2. The Lagrange-mesh method is used for the radial coordinate fl2l] . Since the initial 
state ^'^(O), which decays in time, is excited in the continuum, we add the absorbing 
potentials, —iW{ri) and —iW{r2), to the Hamiltonian of Eq. dTJ). These imaginary 
potentials should be non-zero only where the initial wave function vanishes (ri , r2 > R). 
As the time passes, the two neutrons slowly go away to the infinity (ri , r2 0°), and are 
absorbed by these potentials. 

Calculated El strength distribution, dB(El)/dE, is shown in Fig. [T] The attractive 
neutron-neutron interaction, Vnn, is chosen to be weaker than the one to reproduce the 
n-n scattering length in free space lll3l] . The peak position is around E^ = 0.8 MeV. 
This should be compared to the recent measurement, E^ = 0.6 MeV lUlll . In order 
to reproduce the observed peak position more precisely, we may need to improve the 
Hamiltonian and the ground state of ^^Li. 
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FIGURE 1. Calculated B{Ei) distribution for ' ^Li. 



TDHF: Giant dipole resonances in 




Calculations based on the time-dependent Schrodinger equation of Eq. ^ becomes 
increasingly difficult as the particle number increases. A possible way to avoid this 
difficulty is to resort to the density functional (mean-field) theory. In this section, we 
use the time-dependent Hartree-Fock (TDHF) method with the Skyrme interaction to 
study the dipole response in neutron-rich Oxygen isotope, ^^O. 

Although the TDHF wave packet is capable of describing nuclear dynamics of the 
large-amplitude nature, it is not an easy task to requantize the TDHF trajectories. Thus, 
we focus our discussion on its small- amplitude (perturbative) limit. This is nothing but 
the well-known random-phase approximation (RPA) which has been extensively utilized 
for studies of both high-energy giant resonances and low-energy collective excitations 
|I3 • Here, we calculate it in the time-dependent manner. 

The RPA equation for the linear response with the density-dependent interaction is 
usually derived from the small- amplitude approximation of the time-dependent Hartree- 
Fock (TDHF) theory. Thus, we should be able to obtain the same information from the 
time propagation of a Slater determinant. Let us assume that the state is initially the 
ground (HF) state, then a perturbative external field, Vext(0' is switched on. Each single- 
particle orbital follows the TDHF equation 



where e,- are the HF single-particle energies at the ground state. Decomposition of 
\5\j/i{t)) into normal modes leads to the RPA eigenvalue equations. Instead of doing 
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FIGURE 2. (Top) Calculated time-dependent El moment for O. (Bottom) Calculated photoabsorp- 
tion cross section for ^^O. 



this, we may simply calculate the time-dependent £ 1 dipole moment, 

A A 
i=l i=l 

Then, the Fourier transform of Eq. (fTTT) provides information of eigenenergies and 
transition strength of the RPA normal modes. 

We perform the TDHF calculations using the three-dimensional (3D) Cartesian mesh 
representation. The full Skyrme functional with the SGII parameter set is adopted, which 
preserves the Galilean invariance. The time-dependent method provides an efficient tool 
to calculate the E 1 strength function for a wide range of energy. In fact, to diagonalize 
the RPA matrix in the 3D coordinate space is a time-consuming procedure 1 15. 16,.17ll . 
The small- amplitude TDHF gives a feasible alternative. 

In the top panel of Fig. |2l we show the time evolution of Z)io(0 ^^O' Eq. (fTTll . 
with an instantaneous El external field, VextlO = T?^^io5(0' where rj is an arbitrary 
parameter. The time evolution is calculated up to ? = 20 h/MeV. Here, we do not use the 
absorbing potential. Therefore, the particle continuum is discretized, leading to discrete 
peaks in the strength function. Then, we calculate the Fourier transform of Dio{t) with 
complex energy, E + iT with F = 1 MeV. The imaginary part of the energy gives a 
smoothing width to every peak. In the bottom panel of Fig. |2l photoabsorption cross 
section, calculated in this way, is presented. The main peak position is almost the same 
as that of ^^O Mi Experimental data obtained at GSI suggest a low-energy peak around 
Ex = 9 MeV iTlsIi . The present calculation produces the lowest peak at 13 MeV. This 



discrepancy may be due to the fact that the neutron separation energy is too large in 
the calculation with the SGII parameter set. The main peak position around Ey^ = 20 
MeV has not been confirmed experimentally, however, the result is similar to that of the 
relativistic RPA calculation [19]. 

SUMMARY 

We present the time-dependent method to study nuclear response functions. The con- 
tinuum can be properly treated with the absorbing boundary condition. The method is 
first applied to a three-body model to investigate B{E\) distribution in ^^Li. The calcu- 
lation produces a strong £1 strength at low energy, however the peak position is slightly 
too high compared to the experiment. It seems to be necessary to improve the model 
Hamiltonian in order to reproduce the peak position quantitatively. We next apply the 
time-dependent Hartree-Fock analysis to giant resonances. We show E\ resonances in 
^^O, using the SGII Skyrme functional. The main peak is consistent with the relativistic 
RPA calculation. The position of the low-energy peak is calculated to be higher than 
the experiment by about 4 MeV. The time-dependent method is an efficient tool to in- 
vestigate nuclear response function, especially for its bulk structure in a wide range of 
energy. 
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